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Abstract 

^ We compute the intrinsic volumes of the cone of positive semidefinite matrices over the 

real numbers, over the complex numbers, and over the quaternions, in terms of integrals 
related to Mehta's integral. Several applications for the probabilistic analysis of semidefinite 
programming are given. 
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i — i 1 Introduction 



The classification of symmetric cones, also known as self-scaled cones, i.e., closed convex cones, 
which are self-dual, and whose automorphism group acts transitively on the interior, is a well- 
known result. It says that every symmetric cone is a direct product of the following basic 
families of symmetric cones: 



• the Lorentz cones, which have the form {x G W 1 \ x n > (x\ + . . . + a^-i) 1 / 2 }, 

• the cones of positive semidefinite matrices over the real numbers, the complex numbers, 
or the quaternions, 



in 
o 

^ • the single (exceptional, 27-dimensional) cone of 3 x 3 positive semidefinite matrices over 

the octonions. 

This result follows from the theory of Jordan algebras, which is intimately related to the theory 
of symmetric cones, cf. [15]. 

On the other hand, self-scaled cones form the basis of interior-point methods in convex 
optimization. This has been observed in the mid '90s, cf. [31, 32, 23, 16], cf. also the book [35] 
and the survey article [24]. A detailed understanding of these cones, in particular of its statis- 
tical properties, is thus of fundamental importance. By 'statistical properties' we mean basic 
probabilities like the probability that a random convex program is feasible, etc. In the case of 
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linear programming (LP) these fundamental statistical properties are well understood. For ex- 
ample, the probability that a random LP is feasible has been computed by Wendel in 1962 [43], 
and further elementary probabilities can be deduced from this, cf. [12, Thm. 4]. In fact, the 
probability that a random convex program is feasible can be expressed in terms of the intrinsic 
volumes of the reference cone. The intrinsic volumes of a closed convex cone form a discrete 
probability distribution that to some extent captures the statistics of the cone. In the case 
of linear programming the reference cone is the positive orthant, and its intrinsic volumes are 
given by the symmetric binomial distribution, cf. Section 2.1. As for second-order cone pro- 
gramming, the reference cone is the Lorentz cone, whose intrinsic volumes are also, basically, 
given by the symmetric binomial distribution (see Remark 3.7 for the details). 

We give in this paper, apparently for the first time, an explicit formula for the intrinsic 
volumes of the cone of positive semidefinite matrices over the real numbers, over the complex 
numbers, and over the quaternions. The resulting formulas involve integrals that are related 
to Mehta's integral. It remains to give closed formulas and/or derive asymptotics for these 
integrals, but we see this work as the first step for understanding the intrinsic volumes of the 
symmetric cones, and thus for understanding fundamental statistical properties of semideflnite 
programming. To illustrate the significance of the intrinsic volumes (and of its local versions, 
the curvature measures), we will give several corollaries, which describe some interesting prob- 
abilities about semidefinite programming (SDP) in terms of these integrals. In particular, we 
obtain a closed formula for the probability that the solution of a random SDP has a certain 
rank. To the best of our knowledge, this is so far the first result making concrete statements 
about the above-mentioned probability, a question which is by now at least 15 years old, cf. [4]. 

Another interesting aspect, which deserves further investigation, is the observation that 
there seems to be a connection between the curvature measures of the cone of positive semidef- 
inite matrices and the algebraic degree of semidefinite programming. This degree has been 
defined in [33], cf. also [11], and some remarkable parallels to the curvature measures can be 
established, cf. Remark 3.6. In fact, the authors also describe in [33, Sec. 3] an experiment to 
empirically analyze the rank of the solution of a random semidefinite program. But their (ex- 
perimental) results are not comparable to our formulas, as they choose a different distribution 
to specify a random semidefinite program. 

The organization of the paper is as follows. In Section 1.1 we describe the connection 
between the specific curvature measures, for which we will give explicit formulas in Section 3.3, 
and some fundamental statistics in semidefinite programming. Although we defer the formal 
definition of the curvature measures and the intrinsic volumes to Section 3.1, we give the 
SDP application at this early stage to motivate the subsequent sometimes technical sections. 
Section 2 is devoted to the applications of the kinematic formula in the context of convex 
programming. This section is independent from the rest of the paper, but it shows the usefulness 
of the kinematic formula and further motivates the computation of the intrinsic volumes of the 
cone of semidefinite matrices. In Section 3 we give the main result of this paper, which are 
the formulas for the intrinsic volumes /curvature measures of the cone of positive semidefinite 
matrices. Section 4 is devoted to the proof of the main result. 

Since we need in Section 2 a specific form of the kinematic formula, which is not easy 
to trace in the literature, we will describe in the appendix the additional concept of support 
measures. The most general form of the kinematic formula in the spherical setting is stated in 
terms of these measures. We will describe how the specific kinematic formula that we use can 
be derived from this general result. 
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1.1 Applications in convex programming 

We consider the following forms of convex programming. Let £ be a finite-dimensional Eu- 
clidean space with inner product (.,.): £ x £ — >■ M. Furthermore, let C C £ be a closed convex 
cone, i.e., a closed set that satisfies £ C and Xx + fiy £ C for all x,y £ C and A, /i > 0. The 
classical convex programming problem (with reference cone C) has the inputs a±, . . . , a m , z £ £ 
and &i, . . . , b m £ M, and consists of the task 

maximize (z, x) (CP) 
subject to (ai,x) = bi , i = 1, . . . ,m, 
x £ C, 

which is to be solved in x £ £ . We also define a homogeneous version, which is easier to analyze. 
This has only the inputs ai, . . . , a m , z £ £, and is again to be solved in x £ £ 

maximize (z, x) (hCP) 
subject to (aj, x) = , i = 1, . . . , m, 
x £ C , \\x\\ < 1. 

The space £ is endowed with the (standard) normal distribution M{£). Choosing an or- 
thonormal basis of £ so that we have an isometry (p: £ — > M. d , d = dim£, a random element 
x £ £ is normal distributed iff the components of(p(x) are i.i.d. standard normal, i.e., in jV(0, 1). 
We say that an instance of (CP) is a (normal) random program if the inputs a\, . . . , a m , z are 
i.i.d. in N{£) and the inputs bi,...,b m are i.i.d. in W(0, 1). Analogously, we speak of a random 
instance of (hCP) if the inputs oi, . . . , a m , z are i.i.d. in N{£). 

For the analysis of the problems (CP) and (hCP) we use the following notation: We denote 
the feasible set of (CP) and (hCP) by 

J"(CP) := {x £ C | Vi : {a, h x) = 6J, J"(hCP) := {x £ C \ Vi : (a», x) = 0, ||x|| < 1}. 

By val we denote the value of (CP) or (hCP), 

val(CP) := sup{(z,x) | x £ -F(CP)}, val(hCP) := sup{(z,x) | x £ ^(hCP)}. 

Finally, we denote by Sol the solution set of (CP) or (hCP), 

Sol(CP) := {x £ JF(CP) | (z,x) = val(CP)}, Sol(hCP) := {x £ JT(hCP) | (z,x) = val(hCP)}. 

Note that val(hCP) is in fact a maximum, as the set ^(hCP) is always compact and contains the 
origin. In general, this is not the case for the affine version (CP). The feasible set J-"(CP) may 
be unbounded, and the value val(CP) may be oo in which case we say that (CP) is unbounded. 
Also, the feasible set .F(CP) may be empty, so that val(CP) = sup0 which is — oo by definition. 
In this case we say that (CP) is infeasible. If the solution set Sol(CP) only consists of a single 
element, then we denote this by sol(CP). So writing xq = sol(CP) means that Sol(CP) = {xq}. 
We use a similar convention for Sol(hCP). Well-known results from convex geometry, see for 
example [37, Thm. 2.2.9], imply that almost surely Sol(hCP) and Sol(CP) are either empty or 
consist of single elements. 
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For random instances of (CP) and (hCP) we can express certain statistics in terms of the 
intrinsic volumes and the curvature measures of the reference cone. The intrinsic volumes of 
the closed convex cone C C £ are nonnegative numbers Vb(C), . . . , Vd(C), d = dim£, which 
add up to one, Vq(C) + . . . + Vd(C) = 1. To describe the curvature measures, we make the 
following definition (cf. Section 2.1 for a discussion of this concept). 

Definition 1.1. Let £ be a finite-dimensional euclidean space, and let £§(£) denote the Borel 
a- algebra on £ . Then we call 

®{£) := {M G &{£) | VA > : AM = M} (1) 

the conic (Borel) a-algebra on £. 

It is easily seen that indeed satisfies the axioms of a u-algebra. The curvature measures 
of C are measures $ (C, .), . . . , $ d (C, .) : &{£) -»• M+, which satisfy $j{C,£) = Vj(C). So the 
curvature measures localize the intrinsic volumes. We give an illustrative characterization of 
the intrinsic volumes and the curvature measures of polyhedral cones in Section 2.1, and we 
shall provide the definition for the general case in Section 3.1. In Section A in the appendix 
we will also explain the support measures, which, among other things, will further justify our 
definition oi38{£). 

The following two theorems describe the statistics of (hCP) and (CP). 

Theorem 1.2. We have the following probabilities for random instances of (hCP) 

m d 

Prob [sol(hCP) = 0] =J2Vj(C) , Prob [sol(hCP) € M] = ^ $j(C,M) , (2) 

j=0 j=m+l 

where M £ 3S{£) with M. Furthermore, if C is not a linear subspace, then 

m—l 

Prob [^(hCP) = {0}] =2- W 7 )- ( 3 ) 

3=0 
j=m—l mod 2 

Theorem 1.3. We have the following probabilities for random instances of (CP) 

m— 1 d 

Prob [CP infeasible] = ^ V 3 {C) , Prob [CP unbounded] = ^ Vj(C) . (4) 

j=0 j=m+l 

Furthermore, for M G 3§(£) we have 

Prob [sol(CP) € M] = & m {C, M) , (5) 
and Prob [sol(CP) G M A val(CP) > 0] = Prob [sol(CP) € M A val(CP) < 0] . 

Remark 1.4. The intrinsic volumes of the positive orthant are given by the symmetric 
binomial distribution Vj(W[) = ( d ) /2 d , cf. Remark 2.2. Plug ging in these values in Theorem 1.3 
yields the statistics of linear programming, which are repeatedly computed, cf. [2, 41, 12]. 



4 



We finish this introduction with a discussion of the special case of semidefinite programming. 
For this we need to set up some notation first. Throughout the paper we use the parameter 
(3 G {1, 2, 4} to indicate if we are working over the real numbers R, over the complex numbers C, 
or over the quaternions H. In particular, we denote the ground (skew) field by Fp, i.e., 

Fi := R , F 2 := C , F 4 := H . 

We consider these (skew) fields with the natural identifications 

R C C = R[i] C H = C[j, k] = R[i, j, k] , 

where i 2 = j 2 = k 2 = — 1, and with the well-known quaternion multiplication rule found by 
Hamilton in 1843. Recall that we have the (linear) conjugation map 7: H — > H, given by 1 = 1, 
i = — i, j = — j, k = — k, and its restriction to C. Furthermore, we call 3ft: F^ — > R, z 
the canonical projection on R. 

The set of (n x n)-Hermitian matrices over F^ shall be denoted by 

Her^ := {A G Ff n \ A* = A} , 

where A* = (aji) for A = (aij). This is a real vector space of dimension 

dp, n ■= dimHer^ = n + /3 ■ (™) . (6) 

Throughout this paper we regard Her^ as a euclidean vector space with the inner product 
given by A • B := ^t(tr(A^ B)), where A,Be Her ( g j7l , and tr(A) := an + . . . + a nn denotes the 
trace. The normal distribution in Herg in with respect to this inner product is called the Gaus- 
sian Orthogonal/Unitary /Symplectic Ensemble (GOE/GUE/GSE), according to (3 = 1,2,4, 
respectively. We use the short notation G/3E for this distribution. 

For A 6 Her^ and x G F^ the product Ax lies in R, and an element A 6 Her^ is called 
positive semidefinite iff x^ Ax > for all x G F^. The set of all positive semidefinite elements 
in Herg in is a closed convex cone, the cone of positive semidefinite matrices overFp, which we 
denote by 

Cp, n = {A G Her^ | Vx G F^ : x^ Ax > 0} . (7) 

The cone Cp jTl has a natural decomposition according to the rank of the matrices (cf. [45] for 
the quaternion case) 

n 

Cp,n = |J Wp, n , r , Wp, n , r ■= {A G Cp, n \ rkA = r} . (8) 

r=0 

For the jth curvature measure of Cp :Tl evaluated at the set of its rank r matrices we write 

*j03,n,r) :=®j{Cp, n ,Wp, n , r ) . (9) 
The decomposition (8) of the cone Cp )7l into the rank r strata yields the formula 

n 

^(y^fjIA^r), j = 0,...,dp, n . (10) 

r=0 
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Figure 1: The intrinsic volumes of C 4j 3 and their decompositions in curvature measures. The 
small numbers indicate the contributions of the ranks. The probabilities from Corollary 1.5 for 
m = 6 are also indicated. 



Semidefinite programming is now the following specialization of convex programming 
£ = Rerp, n , C = Cp, n , (x,y)=XmY, M{£) = G/3E , 
and we obtain the following programs 

(SDP^) 



max. Z • X 

s.t. Ai • X = bi 

XhO, 



max. Z • X 
s.t. Ai»X = 

X h 0, ||X|| < 1, 



(hSDP^) 



where X y is the usual notation for X £ Cp n . 

Specializing Theorem 1.3 yields the following corollary for semidefinite programming. 

Corollary 1.5. We have the following probabilities for random instances of (SDP^) 

m— 1 d 

Prob [SDP^ infeasible] = ^ Vj(Cp in ) , Prob [SDP^ unbounded] = ^ Vj(C^ n ) . (11) 

j=0 j=m+l 



Furthermore, for < r < n we have 

Prob [rk(sol(SDP^)) = r] = $ m {P,n,: 



(12) 



See Section 3.3 for explicit formulas for Vj(C/3^ n ) and $>j(/3,n,r). For /3 = 4, n = 3, m = 6 
the probabilities from Corollary 1.5 are shown in Figure 1. 



2 Applications of the kinematic formula 

The goal here is to provide the proofs for Theorem 1.2 and Theorem 1.3. In Section 2.1 we 
first introduce the notion of intrinsic volumes and curvature measures for polyhedral cones; the 
case of general closed convex cones is deferred to Section 3.1. The kinematic formula will be 
presented in Section 2.2. In Section 2.3 and Section 2.4 we prove Theorem 1.2 and Theorem 1.3, 
respectively, by means of the kinematic formula. 
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2.1 Intrinsic volumes of polyhedral cones 

The intrinsic volumes of closed convex cones are usually defined in the spherical setting, which 
is obtained by intersecting the cone with the unit sphere. We will give the definition of the 
spherical intrinsic volumes, which is a bit technical, in Section 3.1. However, we would like to 
mention at this point the close relationship between the conic cr-algebra £$(£)-, that we defined 
in (1), and the Borel algebra on the unit sphere S{£) = {x G £ \ \\x\\ = 1}. Namely, we 
have the decomposition S(£) = & (£) U where Jb(£) := {M G S(£) | G M} and 

£%$(£) := {M G &(£) | G" M}. A moment of thought reveals that the mappings 

J (£) -> &(S(£)) , M ^ M n S{£) , J (£) -»• @{S{£)) , M^Mn 5(f) 

are bijections, i.e., we may identify both ^o(£) and £$%{£) with the Borel algebra 3$(S{£)). 
So it might seem that our definition of SS(£) is superfluous, or overly pedantic. But in fact, 
the use of £${£) is not only convenient, as we will see in the course of this section, but also 
valuable in the context of the support measures, that we will describe in Section A. In the 
following paragraphs we will give an illustrative characterization of the curvature measures and 
the intrinsic volumes for polyhedral cones, i.e., intersections of finitely many closed half-spaces. 

If C C Mr is a closed convex cone, we denote by C := {x G R d \ My G C : (x,y) < 0} the 
dual cone of C in M. d . (Occasionally, we will also use the notation dual(C) := C, if this is more 
convenient.) The most important cones used in convex programming are self-dual, i.e., C = —C; 
it is well-known, cf. for example [8, §11.12], that Cp j7l is self-dual, i.e., dual(Cg )n ) = —Cp^ n . 

A supporting hyperplane H of C is a hyperplane such that C lies in one of the closed 
half-spaces bounded by H. The intersection H n C is called a face 1 of C. 

A polyhedral cone C C R d is the intersection of finitely many closed half-spaces bounded 
by linear hyperplanes. The boundary of the cone C decomposes in the disjoint union of the 
relative interiors of its faces. More precisely, we have 

C = (J f ^F, T := {relint(C n u ± ) | v G C} , 

where v 1 - := {x G M d | (x,v) = 0}. Let Tj := {F G T \ dim(spanF) = j} denote the set of 
(the relative interiors of) the j'-dimensional faces of C, j = 0, 1, . . . , d. 

For a spherical Borel set M s G ^(S^ -1 ) we can write the (d — l)-dimensional normalized 
Hausdorff volume, which we shall denote by rvol(M s ), in the form 

rvol(ikP) = Prob \x G M] , (13) 

where Af(0,ld) denotes the standard normal distribution on R d , and M := {Ax | A > 0, x G 
M s } G J(R d ). Clearly, the distribution Af(0,I d ) may be replaced by any other orthogonal 
invariant distribution [i on M. d , which satisfies //({0}) = 0. 

Denoting by lie: R d — > C, x i— > argmin{||x — y\\ \ y G C} the canonical projection on C, 
the intrinsic volumes of C are given by 

V ^ = E ^° n b r , i U c(x) €F], j = 0,l,...,d. (14) 



1 Some authors differentiate between faces and exposed faces, cf. for example [37]. We do not make this 
distinction, as for those cones in which we are interested both notions coincide. 
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Note that V d (C) = rvol(C n S^ 1 ) and V (C) = rvol(C n S^ 1 ). 

For M G the curvature measures of C evaluated in M are given by 



^■(C,M)=V Prob [n c (i)GFnM], j = 0,l,...,d. (15) 
Note that for j G {0, d} we have 

*,(C,M)^vo,<CnMnS->, ...a,^*'-^ 5 '- 1 ' Hl M M (16) 

Additionally, we define Vj(C) := and ^(C, M) := for j > d. 

One could use the formulas (14) and (15) to define the intrinsic volumes and curvature 
measures for general closed convex cones, using an approximation procedure. But a more useful 
definition is via a spherical version of Steiner's formula for the volume of the tube around a 
convex set. We will describe this in Section 3.1. 

The following well-known facts about the intrinsic volumes and the curvature measures may 
be verified easily for polyhedral cones using the above characterizations of Vj and &j in (14) 
and (15). 

Proposition 2.1. 1. Interpreting C C ]R d as a cone in M. d ' with d! > d does not change the 
intrinsic volumes nor the curvature measures. 

2. The intrinsic volumes and the curvature measures are nonnegative, and ~}2 d j = QVj{C) = 1 
i/CC R d . We have Vj{W) = <% the Kronecker delta. 

3. The curvature measure &j(C, .) is concentrated on C, that is, §j(C,M) = $j(C,MnC). 
Furthermore, we have &j(C,C) = Vj{C). 

4- The intrinsic volumes and the curvature measures are invariant under orthogonal trans- 
formations, i.e., for Q G 0(d) we have Vj(QC) = Vj{C) and $j(QC,QM) = $j(C,M). 

5. For the intrinsic volumes of the dual cone we have Vj(C) = V d -j(C). 

6. If C\,C2 are closed convex cones, then Vj(C\ x C2) = Si=o^»(^i) ' ^j-ii^)- I n other 
words, the intrinsic volumes of a product arise as the convolution of the intrinsic volumes 
of the components. 



7. If W C M. d is a linear subspace of codimension m and Tlw '■ K d — > W the orthogonal 

-j+m ( 



projection onto W, then^j(U w (C),U w (M)) = <S> j+m (C+W ± , M+W^) for M G 3g(R d ). 



8. The probability that the projection of a Gaussian vector lies in M 6 38 (Br) is given by 

Prob [U c (x) G M] = E -=0 ^(C, M). □ 
x€N(0,I d ) 

Note that the self-duality of Cp jTl implies 

V j (Cp, n ) = V dptn . j (C fi , n ). (17) 

An important but nontrivial property of the intrinsic volumes is the following consequence 
of the Gauss-Bonnet Theorem: For a closed convex cone C C R d 

v 1 (C) + v 3 (C) + v 5 (C) + ... = l-x(cn s d - 1 ) , (is) 
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where \ denotes the Euler characteristic, cf. [20, Sec. 4.3] or [38, Thm. 6.5.5]. Note that 
x(C fl S^ 1 ) = 1 if C is a closed convex cone which is not a linear subspace. In this case we 
have 

£ Vj(C) = Y, Vj(C) = \ . (19) 

j even j odd 

Remark 2.2. An important example for a polyhedral cone is the positive orthant M.+ . Its 
intrinsic volumes are easily computed using the product rule (6) in Proposition 2.1: We have 
Vo(M+) = Vi(M+) = \, i.e., the intrinsic volumes of a 1-dimensional ray form a symmetric 
Bernoulli distribution. Hence, the intrinsic volumes of the positive orthant M.+ = R + x . . . x 
M + are the ci-times convolution of the Bernoulli distribution, i.e., the symmetric binomial 
distribution Vj(R%) = ( d j )/2 d . 

2.2 The kinematic formula 

Euclidean versions of the kinematic formulas are well-known, cf. for example the survey ar- 
ticle [28] and the references given therein. Spherical versions, on the other hand, are less 
well-known. The formulas we need in this paper may be derived from a general version due 
to Glasauer [20], cf. also [21] or [28, §2.4] for (short) summaries. As the literature for spher- 
ical intrinsic volumes is sparse and known results sometimes hard to find, we will describe in 
the appendix how the kinematic formulas, as we state them here, are derived from Glasauer's 
result. 

The uniform probability distribution on the set of fc-dimensional subspaces of M. d , the Grass- 
mann manifold Gr^fc = {W C R d | W linear subspace of dimension k}, is characterized as 
the unique probability distribution, which is invariant under the action of the orthogonal 
group 0(d). Loosely speaking, a dimensional subspace is drawn uniformly at random, if 
every subspace 'has the same probability'. This distribution is for example obtained as the 
push-forward of Gaussian matrices (of appropriate format) via taking the kernel or the image. 

Theorem 2.3 (Kinematic formula). Let C C W 1 be a closed convex cone, and let W C R rf 
be a uniformly random subspace of codimension m G {1, . . . , d — 1}, i.e., d\m.W = d — m. 
Furthermore, let M 6 ^(M d ) be such that M C C ' , and let Hw denote the orthogonal projection 
on W. Then we have for the random intersection C CiW 

E[$j(CnW,MnW)] =$ m+j (C,M) , for j = l,2,...,d-m , (20) 

E[V (CnW)] =V (C) + V l (C) + ... + V m (C) , (21) 

and for the random projection Uyy(C) we have 

E[Q j (IL w (C),n w (M))]=* j {C,M), forj = 0,l,...,d-m-l, (22) 
E[V d . m {Ii w {C))] = Vd_ m (C) + V d - m+1 (C) + ... + V d (C) . (23) 

Using the expression of the Euler characteristic given in (18), one obtains from (20) a 
corollary about the probability that the random intersection C n W is the zero set {0}. 

Corollary 2.4. Let C cR d be a closed convex cone, which is not a linear subspace. Then for 
W C W 1 a uniformly random subspace of codimension m 

Prob [CP W = {0}] = 2 • (y m _i(C) + V m - 3 (C) + V m - 5 (C) + ...). (24) 
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Proof. The Euler characteristic x^nVFflS"^ 1 ) vanishes if CnW = {0} and equals 1 otherwise, 
provided C n W is not a linear subspace. Moreover, the intersection C PiW is almost surely 
not a linear subspace. Therefore, we may write the probability for the event C n W / {0} as 
an expectation and apply the kinematic formula 

Prob [c n w ± {o}] = E[ x (Cnwns d - 1 )] ( =2- E[^-(crw)] ( =2- y m+j (C) . 

j odd j odd 

As the intrinsic volumes with even/odd indices add up to \ by (19), we obtain 

Prob [C n W = {0}] = 1 - 2 ■ V m+j(C) ( = 2 • £ V m ^(C) . □ 

j odd j odd 

2.3 Statistical properties of (hCP) 

We introduce the following notation. A subset if C S d ~ l of the (d — l)th unit sphere is called 
spherically convex iff it is of the form if = C n S d '\ where C C R d is a closed convex cone. 
We define the dual set of if via K := C D S^ -1 . Furthermore, we denote the projection map 
onto if by 

s^if ^if , n^(p) := linear 1 -n c (p) , (25) 

where lie denotes the canonical projection onto C. The angle d(p, q) = arccos(p, g) between 
p,q G S^ -1 defines a metric on S^ 1 . We use the notation d(p,K) := min{d(p,q) | g G X}. 
Note that the dual set K may be characterized as if = {p G S"^ 1 | ci(p, if) > §}• 

To simplify the notation we adopt the following convention: If we maximize a function / 
over a set M, then we denote Argmax{/(x) | x G M} := {x G M | f(x) = sup{/(y) | y G M}}. 
If the set Argmax{/(x) | x G M} consists of a single element only, then we denote this element 
by argmax{/(x) | x G M}. Similarly for Argmin and argmin. Note that for v £ S^ 1 \ K we 
have 

Hk{v ) = argmin{<i(p, u) | p G If } = argmax{(p, v) \ p G if} . (26) 

We will see that the homogeneous case (hCP) is easily reformulated in such a way that the 
kinematic formula yields the proof of Theorem 1.2. The key observation is made in the next 
lemma, cf. Figure 2. 

Lemma 2.5. Let C C~R d be a closed convex cone, K :=CC\S d ' 1 , and let B C R d denote the 
closed unit ball. Then for v G S^ 1 

argmax{(,,x)|,GCn5} = (^ W (27) 

10 p£ mt(if ) . 



Proof. For v G" K there exists p G if such that d(p,v) < |, i.e., (p, u) > 0. It follows that 

Argmax{(w, x) \ x G C fl 5} = Argmax{(i>,p) | p G if} ^ {IIr-(v)} . 

On the other hand, if u G int(if ), then d(p, v) > | for all p G if , i.e., (p, u) < for all p G if . 
It follows that in this case Argmax{(f , x) | x G C fl 5} = {0}. □ 
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(a) veK 



(b) v&{K\JK) 



(c) v G K 



Figure 2: An illustration of Lemma 2.5. 



The problem (hCP) can now be phrased in the following form: We have the closed convex 
cone C in d-dimensional euclidean space £. This cone is intersected with the closed unit ball 
B{£) := {x G £ | ||x|| < 1} and with the linear subspace W := {x G £ \ (a\,x) = ... = 
(a m , x) = 0}. In other words, we have 

J-(hCP) =Cn Wf)B(£) . 

If the Oj are from the normal distribution N(£) then W has almost surely codimension m, 
and is uniformly distributed among all (d — m)-dimensional subspaces of £. So we may 
assume w.l.o.g. that W is a uniformly random (d — m)-dimensional subspace of £ . 
We are now in a position to apply the kinematic formula. 

Proof of Theorem 1.2. The linear functional z to be minimized in (hCP) may be replaced by 
its orthogonal projection z on W, as this does not change the value of the functional on W. For 
fixed W we thus obtain a conditional distribution for z, which, by the well-known properties of 
the normal distribution, is again the normal distribution (on W). Hence the probability that 
the origin is the solution of (hCP) is given in the following way 

Prob [sol(hCP) = 0] = Prob Probr a rgmax{(z,x) | x G W n C D B(£)} = 0] 

a\,...,a m ,z ' W z 

m 

{ = ] Prob Prpb[z G dual(W n C)] ( = E [V (W n C)] ( = ^ Vj(C) , 
w w i=o 

which shows the first claim in (2). As for the second claim in (2), let Hc w denote the projection 
onto C w := C n W. Then we obtain for M G $(£) such that M 

Prob [sol(hCP) G Ml = Prob Prob[argmax{(z, x) I x G If flCnM)} G Ml 

ai,...,a m ,z L J iy 2 L ' J 



(27) 

For fixed W we obtain from Proposition 2.1(8) 



Prob Prpb[n CH/ (z) G M] . 



d—m 

Frpb[n Cw (z) G M] = ^j(cnw,M) . 

3=1 

For random W we may apply the kinematic formula and obtain 

d—m d—m 

Prob [sol(hCP) G M] = V E [^(CflW.M)] ( = } V $ j+m (C,M) , 

All — )^TH)2 W 
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which shows the second claim in (2). 

Finally, if the cone C is not a linear subspace, then by Corollary 2.4: 

Prob [^(hCP) = {0}] = Prob \C n W = {0}1 

ai,...,a m W 

( = } 2 • (V m -i(C) + V m - 3 (C) + V m - 5 (C) + ...), 
which is (3) and thus finishes the proof of Theorem 1.2. □ 

2.4 Statistical properties of (CP) 

The geometric interpretation of (CP), which is suitable for applications of the kinematic for- 
mula, is slightly more complicated than in the homogeneous case (hCP). The key observation 
is in the following lemma, which reduces the d-dimensional to the 2-dimensional case. 

Lemma 2.6. Let v, w G be such that (v, w) = 0, and let L := spanju , w} denote the plane 
spanned by v and w. Furthermore, let II^: M. d — > L denote the orthogonal projection onto L. 
Then for a closed convex cone C C ~R d we have 

sup{(w, x) I x G C , (w, x) = 1} = sup{(f , x) | x G ITl(C) , (u;, x) = 1} , 
Argmax{(v, x) | x G C , (to, x) = 1} = C D II^ 1 (Argmax{(f , x) | x G n^(C) , (w, x) = 1}) . 

Proof. Let x G M rf be decomposed in x = xi + X2 with xi G L and X2 G L -1 -, i.e., Xi = n^(x). 
Then we have (v,x) = (v,x±) + {v,X2) = (f,xi), and similarly (w,x) = (w,xi). This implies 
sup{(w,x) | x G C, (w,x) = 1} = sup{(i;,xi) | x\ G Hl(C) , {w,x±) = 1}. Analogously, we 
obtain the second claim. □ 

We now discuss the 2-dimensional case. Let v , w G S 1 with (v, w) = 0, i.e., the matrix with 
columns v, w lies in 0(2). The orthogonal group 0(2) is isometric to the disjoint union S 1 U S 1 

7/7 11 1 ) 111 

according to the two possible orientations V and V . To make this explicit we 
define the map 

^. . * r „ -, ^ . . |(f,^) if the orientation is w N/ 1 w , . , 

(^:0(2)^ ( S 1 x{£,r}, Q = (v,w) m- ' ' (28) 

[(f,r) if the orientation is N/* . 

In the following let C C I 2 be a fixed closed convex cone, which is not a linear subspace, 
i.e., C is a wedge with an angle between and it. Figure 3 illustrates where the intersection 
of C with the affine line {x G M 2 | (x, w) = 1} achieves its maximum with respect to the linear 

functional defined by v, where (v,w) G 0(2) according to the orientation w \S v . 

Let the two rays forming the boundary of C be denoted by R\ and i?2- Furthermore, 
depending on v, w, write 

T := {x G C | (w,x) = 1}, val := sup{(f,x) | x G T}, Sol := Argmax{(v, x) | x G F}. (29) 

Assuming that v, w are random vectors, such that (v,w) G 0(2) uniformly at random, it is 
easily seen that only four cases appear with positive probability: The intersection T may be 
empty, the functional v may be unbounded on T, or the solution set Sol consists of a single 



12 



Figure 3: The 2- dimensional situation in the inhomogeneous case. 



\ C / 

dual(C) 

(a) Notation 




S0I.G Ri 



(b) Decomposition of 0(2) = S 1 x {£, r} according to (29). 




Figure 4: Illustration of Lemma 2.7. 



point, which either lies in R± or in R2. In the latter case we again adopt the convention 
to denote the single point by sol. This case distinction corresponds to the decomposition of 
0(2) = S 1 x {£, r} indicated in Figure 4, where we use the isometry cp defined in (28). 
We formulate the following lemma, which is checked easily, cf. Figure 4. 

Lemma 2.7. Let C C M 2 be a closed convex cone, which is not a linear subspace, and let R\ and 
i?2 denote the two rays forming the boundary of C. Then, for uniformly random (v, w) G 0(2), 
we have 



Prob [F = 0] = V (C) , Prob 



val 



00 



where we use the notation from (29) . Furthermore, for M G 



= V 2 (C) , 
l 2 ), we have 



(30) 



Prob 



sol G M 



= $i(C,M), Prob 



sol G M and val > 



= l-^(C,M). □ (31) 



As in the homogeneous case, we will now bring the problem (CP) into a geometric form 
where we can apply the kinematic formula. This requires a few more steps than for (hCP). 

In addition to W := {x G £ \ (ai,x) = ... = (a m ,x) = 0} we introduce the following 
notation 

Wag := {x G £ I (ai,x) = h, . . . , {a m ,x) = b m } , W := span(W afr ) . 

For normal distributed a±, . . . , a m , 61, . . . , b m , the set W^s is almost surely an affine space of 
codimension m and its linear hull IF is a linear space of codimension m — 1. Moreover, the 
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Figure 5: An illustration of the geometric situation. 

probability distribution of the a« and 6j induce for W a uniform distribution on the Grassmann 
manifold of subspaces of £ with codimension m — 1. 

The affine space W a ff has almost surely a positive height min{||x|| | x G W a ff}. We define 
the normalization W° ff of W a ff, which has unit height, via 

Wtf = h' 1 ■ WW , /i = min{||x|| | x G W aff } . 

Additionally, we define the point w G W° ff by the property W° s PI S'(W) = {w}. So we have 
W = W n w;- 1 - and W° ff = it; + W. See Figure 5 for an illustration of these definitions. 

As we are not interested in the specific value of (CP) (provided it is < oo) but only where 
the maximum is attained, we may consider W° ff instead of W a ff , i.e., instead of (CP) we consider 

maximize z ■ x s.t. x G C fl W , (u;, x) = 1 . (32) 

It is easily seen that for fixed W the induced distribution of w is the uniform distribution 
on S{W). 

Without loss of generality, we may replace the functional z by its orthogonal projection z 
on W. For fixed W the induced distribution of z is the normal distribution on W. As z is 
almost surely nonzero, we may define the normalization v := • z G S(W). Finally, we 

denote the plane spanned by v,w by L := span{v, w}. 

We can generate the distribution of (W, L, v, w) induced by the standard normal distributed 
ai, . . . , a m , &i, . . . , b m in the following way: 

1. choose a uniformly random subspace W of £ of codimension m — 1, 

2. choose a plane L QW uniformly at random, 

3. choose v G S(L) uniformly at random, 

4. choose w as one of the points in S(L) n v -1 , each with probability \. 

We may now proceed to the proof of Theorem 1.3. 

Proof of Theorem 1.3. Lemma 2.6 tells us that instead of (32) we may consider the following 
problem in the 2-dimensional plane L 

maximize (v,x) s.t. xGn^(CnW), (w,x) = l. (33) 
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More precisely, using the notation of (29) 

C :=U L (Cr\W) , F:= {x G C | (w,x) = 1} , 

val := sup{(v,x) | x G F} , Sol := Argmax{(v, x) \ x 6 .F}, 

we obtain from Lemma 2.6 that (CP) is infeasible iff F = 0, (CP) is unbounded iff val = 
and 

Soi(CP) = cnn-^Sol) . 



Prob [CP infeasible] = Prob Prob \F = 0] ( = } E \V (U L (C D W))] 

ai,...,am L ■ w r v,w J W L 



We thus obtain by Lemma 2.7 

,:■!()) 

\Kyf mieasiDiej = jttod rroDiy- = v\ 

, * ' W,L v > w 

bi,...,b m 

Applying the kinematic formula twice yields (recall that the codimension of W is m — 1) 

e [v (u L (cnw))] { = ] E[v {cnw)] { = ] v Q (C) + v l (C) + ... + v m ^(C) , 

W,L W 
which proves the first claim in (4). Analogously, we obtain 

Prob [CP unbounded] = Prob Prob [val = oo] ( = } E \V 2 (U L (C n W))] 

«i,...,am,* L w,L v,w 1 J W,L l 1 

bl,...,bm 

{ = ] e [v 2 (c nw) + v 3 (c nw) + ... + v d - m+1 (c n w)] 
w 

{ = ] V m+1 (C) + V m+2 (C) + ... + V d (C) , 

which proves the second claim in (4). 

As for the claim (5), we have for M G (cf. (34) above) 

(sol(CP) G M n W and val(CP) > 0) (sol G n L (M n W) and v^d > 0) . 

This yields 

Prob [sol(CP) G M and val(CP) > 0] = Prob Probfcu" G U L (M n W") and v^I > 

a lr-i a mi2 " £ 1),UI L 

bi,...,6 m ' 

( = } e [i-$ 1 (n L (cn#),n L (Mnw))l . 

Applying the kinematic formula twice finally yields (recall codiml^ = m — 1) 

e [§ • $ 1 (n L (cniy),n L (Mn vf))] ( i } ± ■E[$i(cnf,Mnw)] ( = \ ■ $ m (c,M) 

W,L W 
Analogous arguments yield the claim without the positivity assumption. 
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3 Main result 



In this section we will formulate our main result, which are closed formulas for the curvature 
measures of the cone of positive semidefinite matrices over R/C/H evaluated at the set of rank r 
matrices. We first give the definition of intrinsic volumes and curvature measures of arbitrary 
closed convex cones in Section 3.1. Section 3.2 is devoted to Mehta's and related integrals. 
The notation we establish here will allow us to formulate our main result in Section 3.3. In 
Section 3.4 we will give some examples in small dimensions. 

3.1 Spherical intrinsic volumes 

Central to the definition of the intrinsic volumes and the curvature measures is the notion of the 
(local) tube around a spherically convex set K C S d_1 . For a spherical Borel set M s G ^(S^ 1 ) 
we define the (local) tube around K (in M s ) of radius a G [0, 7r/2) via 

T(K, a) := {p G S^ 1 \ d(p, K) < a] , 
T(K, a; M s ) := {p G T(K, a) | U K (p) G M s } , 

where Hk denotes the spherical projection map, cf. (25). We oppress the dependence on S^ 1 
to keep the notation simple. Recall that we denote the (d— l)-dimensional normalized Hausdorff 
volume on S"* -1 by rvol, cf. (13). 

The following proposition forms the basis for the general definition of the curvature measures 
and the intrinsic volumes. For a proof see for example [26, 5, 36, 29, 20]. 

Proposition 3.1. Let K C S" i_1 be a spherically convex set. Then there exist nonnegative 
measures $o(K, .),<f>i(K, .),... ^^(K, .): m{S d ~ x ) R + such that for < a < vr/2 and 
M s G SB{S d ~ x ) 

d-i 

ivo\T{K,a;M s ) = $ d _i(K, M s ) + ^ $j-i(K, M s ) ■ rvolT^ 1 , a) . (35) 

3=1 

Furthermore, if K±, K2, ■ . . C S 11 " 1 is a sequence of spherically convex sets, which converges 
to K in the Hausdorff metric, then lim^oo Qj-^K^M 8 ) -»• ^_i(K,M s ) for all M s G 

Definition 3.2. Let C C R d be a closed convex cone, and let K = CnS d 1 be the corresponding 
spherically convex set. For 1 < j < d, the function Qj_i(K, .) from Proposition 3.1 is called 
the (j — l)th curvature measure of K. The intrinsic volumes of K are defined by 

Vj-iiK) :=$> 1 - 1 (K,S d - 1 ) . 

Let the conic a- algebra 3§(R d ) be defined as in Section 2.1. For 1 < j < d the jth curvature 
measure of C is the functional &j(C, .) : £$(R d ) — > R+, which is given by 

$j(c,M) -.= ^ j . 1 (cns d - 1 ,Mns d - 1 ) . 

The curvature measure $o(C, •) : 3${R d ) ->• R+ is defined by the scaled Dirac measure 



$o(C,M) := 




^CnS*- 1 ) ifOGM, 
if M . 
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The intrinsic volumes of C are defined by 



Vj{C):=*j{C,C), j = 0,...,d. 

Remark 3.3. Note that we introduce a shift in the index of the curvature measures/intrinsic 
volumes when passing between the cone and its intersection with the unit sphere. We do 
this for several reasons. First of all, the spherical notation is established this way at several 
places [20, 21, 24, 6]. Second, if interpreted correctly, Vj(C) and Vj-i(K) may indeed be seen 
as volumes of sets with "dimension" j and j — 1, respectively. And last but not least, we have 
made the experience that the formulas get nicer when working with the shifted index for the 
cones, cp. the kinematic formulas in Section 2.2 and in Section A. 3 in the appendix. 



3.2 Mehta's and related integrals 

For z = (zi, . . . ,z n ) we denote the Vandermonde determinant by A(z) := ni<i<j<n( z « — z j)- 
For < r < n let x := (zi, . . . ,z r ) and y := (z r+ i, . . . , z n ), so that z = (x, y). We have the 
decomposition 

r n—r 

A(zf = A(xf ■ A(yf • J] J] ( Xi - y 3 f . (36) 

i=l j=l 

We regard the rightmost factor in (36) as a polynomial in x and decompose it into its homoge- 
neous parts. For convenience, we change the sign, and define 



f/3,k(x; y) := ( the x-homogeneous part of ]^[(xi + yj)^ of degree k 
^ i=i j=i 

We can write //3,fc(x; y) in & n explicit form if we rearrange 



n—r 



i=i j=i i=i j=i j=i 



Denoting by the kth. elementary symmetric function, we obtain 

n—r 

ffiAv y) = ((x y- 1 )^) • n yf - (37) 

where (x ® y _1 ) x/3 = (x <g> . • ■ , x <8> y _1 ) , and 

V v ' 

/3-times 

x «g, y -l := ( ^1,...,^,^,...,^,..., _^ j ...,^ > ) GK r(n-r) (3g) 
\yi 2/1 2/2 2/2 J/ra-r 2/n-r/ 

We thus have a decomposition of the /3th power of the Vandermonde determinant into 

/3r(n— r) 

A{zf = A{xf ■ A{yf ■ Y, -i/)> ( 39 ) 

fc=0 

where 2 = (x, y). 
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Definition 3.4. We define for < r < n and < k < (3r(n — r) the integrals 

^(„,r,fc):=^72- / e"^ ■ \A(x)f ■ \A(y)f ■ fa k (x;y) dz , 



(40) 



where z = with i £ R r j £ 1" r , and M™ denotes the positive orthant in W 1 . We set 

Jg(n, r, k) := 0, if A; < or k > (3r(n — r). 

Note that for n = 1 we have 

Jp(l, 0,0) = ^(1,1,0) = I . (41) 
Note also that exchanging the roles of x and y yields the following symmetry relation 

Js(n, r, k) = Jg(n, n — r, j3r(n — r) — k) . (42) 



For r G {0, n} and = we obtain the integrand e 2 • |A(z)|^, which also appears in 
Mehta's integral 



n 



n 



i m 



(43) 



(cf. [18] and the references therein). 

It is well-known that the distribution of the joint probability density function for the eigen- 
values of matrices from the G/3E is given by (cf. [18] and the references therein) 

1 ■ \A(z)f . (44) 



(2tt)»/2 • F n (p/2) 
Using this, we may write the J-integral as an expectation 

Jf,(n,r,k) = F r (P/2).F n - r (P/2). E [l + (A) ■ 1+(B) ■ f p , k {A; B)] , (45) 

BeG/3E(n-r) 

where 1+(A) := 1 if A y and 1+(A) := if A ^ 0, and fp,k{A; B) is the value of /g jfc in the 
eigenvalues of A and i?. 
If we denote 

F+(/3/2):=J^(n,0,0) = ^-^- | e"^ ■ |A(z)|" dz , 



and for n 


= 2, 3 and /3 = 


1,2,4 








1 


2 


4 




2 




1 1 

2 TT 


7T 




3 


1 / 3 3V2\ 
v^\4 2tt y 


3 9 
2 2tt 


540 - 1692 

7T 



(46) 



In fact, the values for i ? ^"(/3/2) are easily computed with any computer algebra system; we ob- 
tained the values for F 3 + (/3/2) differently (cf. Section 3.4 for the details). As for the asymptotics, 
it is shown in [13] that for n — > 00 

F+(/3/2) = 6 (exp (-n 2 • ^ 
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3.3 Formulas for the curvature measures of the rank r stratum 

Recall the values $>j(/3,n,r) of the jth curvature measures of Cp :Tl evaluated at the set of its 
rank r matrices, cf. (9). The following theorem is a main result of this paper. 

Theorem 3.5. Let (3 G {1,2,4}, n 6 N, and < r < n. The curvature measures of the 
semidefinite cone Cp jTl evaluated at the set of rank r matrices, cf. (7)-(9) ; are given by 

^■-K;)'^S^ (47) 

where < j < dp :7l , and dp >n , Jp(n,r,k), F n ((3/2) are defined in (6), (42), and (43), respec- 
tively. 

Note that the intrinsic volumes of the cone Cp :Jl , are given by Vj(Cp jTl ) = Y2r=o ®j(P-> n i r )> 
cf. (10). Using the expression of Jp(n, r, k) from (45), we may alternatively write the curvature 
measures in the form 

lU(A)-l + (B)-U, i - dtA A;B )] . (48) 

BeG/3E(n-r) 

Remark 3.6. 1. The integral Jp(n,r,k) is nonzero iff < k < (3r{n — r). This implies 

n, r) > ^ d^, r < j < d^r + (3r(n - r) . (49) 

One can deduce the implication "=>•" also from Corollary 1.5, which provides a different 
characterization of $j(/3, n, r). Namely, it is well-known (at least in the real case (3 = 1), 
that a random instance of (SDP^) almost surely satisfies dp jT < m < dp^ r + (3r(n — r), 
where r denotes the rank of the solution of (SDP^). These inequalities are known as 
Pataki's inequalities, cf. [4, 34, 33]. 

2. Note that the relation Jp(n,r,k) = Jp(n,n — r,(3r(n — r) — k), cf. (42), implies the 
symmetry 

&j(P,n,r) = <S> d p n -j(f3,n,n-r) . (50) 

This is a refinement of the duality relation Vj(Cp tn ) = Vd p n -j(Cp >n ), which follows from 
the self-duality of Cp jTl , cf. Section 2.1. 

3. Both of the above properties of &j(/3, n, r) also hold for the algebraic degree of semidefinite 
programming, cf. [33, Prop. 9]. We conjecture a deeper reason for this coincidence, which 
should be interesting to explore further. 

Remark 3.7. With Theorem 3.5 we have formulas for the intrinsic volumes of almost all 
symmetric cones: Recall that the characterization theorem of symmetric cones says that every 
symmetric cone is a direct product of Lorentz cones C n := {x G M n | x n > (x? + ... + 4_ 1 ) 1 / 2 }, 
the cones Cp >n , and the exceptional 27-dimensional cone of positive semidefinite (3 x 3)-matrices 
over the octonions, cf. [15]. The intrinsic volumes of C n can be shown to be (cf. for example [7, 
Ex. 2.15]) 

f(n-2)/2\ 



VAC n ) = " U ~ 1] ,i , for 1< j < n - 1 
j v / 2 n / — j — 



/(n-2)/2\ 

V (n = V n (C n ) = K -^^- ■ 2 F 1 (1, \; ^; -1) 



f{n-2)/2\ 
v -1/2 I 

2"/ 2 



for n — y oo 
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where = r ( y+ ±-j^x-y+i) ' anc ^ 2 -Pi denotes the ordinary hypergeometric function (cf. [1, 
Ch. 15]). Furthermore, we have for direct products the simple convolution rule stated in 
Proposition 2.1(6). With Theorem 3.5 we have thus formulas for the intrinsic volumes of 
all symmetric cones, which do not have the exceptional 27-dimensional cone as one of its 
components. 

It is interesting to see how the formulas for the curvature measures in Theorem 3.5 fit into the 
well-known framework of random matrices from the Gaussian Orthogonal/Unitary/Symplectic 
Ensemble. In the following remark we provide such a connection, which may be interpreted as 
a "sanity check" of the formula for <&j(/3,n, r) given in the main Theorem 3.5. 

Remark 3.8. It is easily seen that the projection map 11^ n : Her^ jn — > Cp )TL simply replaces 
the negative eigenvalues of a matrix A £ Her^ n by 0. It follows that a full rank matrix 
A G Herg n has exactly r positive eigenvalues iff the projection Ylc /3n (A) lies in Wp^ r , cf. (8). 
Using Proposition 2.1(8), we obtain that the probability that a matrix A from the G/3E has 
exactly r positive eigenvalues is given by 

Prob [A has exactly r positive eigenvalues] = ^ n , r ) ■ (51) 

3=0 

On the other hand, using the formula (44) for the distribution of the joint probability density 
function for the eigenvalues of matrices from the G/3E, we obtain 

1 fn\ f IHI 2 

Prob \A has exactly r pos. eigenvalues! = - — : — — — • ( ) • / e ~ ■ \A(z)r dz, 

1 J (2tt)"/ 2 • F n (p/2) \r) J 

R r + xR"T r 

(52) 

where := — M.+. And indeed, (52) coincides with (51), as is seen by the following computation 



(52) (3 = 9) 



3=0 



(2tt)»/2 • F„(/3/2 



f3r(n— r) 



Pr(n-r) . N , dp, n / \ T , . , x rf /3,n 

(40) ^ n\ Jp(n,r,j) ^ / n\ Jp{n,r, j - dp/ (47) ^ 

= g g • - g g ■ - £w> ■ 

3.4 Examples in small dimensions 

In this section we give the values of the curvature measures &j(/3, n, r) and the intrinsic volumes 
of Cp :n for dimension n = 1, 2, 3. 

The case n = 1 is of course trivial, and we only mention it for the sake of completeness: 
we have & (p, 1,0) = V (Cp/ = $i(/3, 1, 1) = Vi{Cp/ = \, cf. (41). The remaining values for 
n = 1 are zero. 
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Table 1: The values of $j(/3,3, 1). 



As for the nontrivial cases n = 2, 3, recall that we have the symmetry relations $j(/3, n, r) = 
n _j(/3, n,n — r) and Vj(Cp iTl ) = n —j(Cp jTl ), cf. (50). This halves the number of values we 

have to determine. For n = 2 the J-integrals are easily computed with any computer algebra 

system. We obtain for the curvature measures 

*o(l,2,0) = !-^, $i(l,2,l) = ^ 

$o(2,2,0) = \- ± , #i(2,2,l) = | , $2(2,2,1) = \ 

$o(4, 2, 0) = \ - & , $i(4, 2, 1) = | , $ 2 (4, 2, 1) = £ , $ 3 (4, 2, 1) = \ 

(the remaining values for n = 2 are either 0, cf. (49), or they are obtained from the above values 
via the symmetry (50)). 

As for the case n = 3, the J-integrals Jg(3, r, fc) for r 6 {1, 2} are easily computed with any 
computer algebra system. We then obtain the value Jg(3, 0, 0) = Jp(3, 3, 0) by first computing 
Vi(Cp t 3), . . . , 3 _i(C / g i 3), which only depend on Jp(3, r, k) for r G {1, 2}. By the self-duality 
of Cp >n we have Vo(Cp t 3) = V^^ 3 (Cp^). Combining this with Proposition 2.1(2) we obtain 
2^0(^,3) = 1 - (Vi(Cp\ 3 ) - V 2 (Cp, 3 ) - ... - ^^3-1(^,3)). From Theorem 3.5 we thus obtain 

^(3,0,0) = F 3 (/3/2)-y (C /3 ,3) 

= F 3 (/3/2) ■ \ ■ (1 - Vi(C A 3) - F 2 (C W ) - ... - ^,3-1(^,3)) • 

The resulting values of the curvature measures are 

$o(l,3,0) = \ - § , $o(2,3,0) = $o(4,3,0) = 1-^, 

and the values $j(/3, 3, 1) are given in Table 1. (The remaining values for n = 3 are obtained 
via the symmetry (50)). 

The values of the intrinsic volumes of Cp >n , n = 1,2, 3, are summarized in Table 2. 

4 Proof of the main result 

In this section we provide the proof of the main Theorem 3.5. We first describe in Section 4.1 
how the intrinsic volumes and the curvature measures may be expressed in terms of curvature 
by stating Weyl's classical tube formula [44] and a generalization, which holds for a larger class 
of cones. In Section 4.2 we state some general facts about the orthogonal/unitary/ (compact) 
symplectic group, that we will use in Section 4.3 for the proof of Theorem 3.5. 
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Table 2: Intrinsic volumes of Cp >n for n = 1,2,3 (the missing entries for £2,3 and 6^3 are 
obtained via Vj(Cp in ) = V df3n -j(Cp, n )). 

4.1 Expressing intrinsic volumes in terms of curvature 

From the characterizations (14) and (15) one easily obtains elementary formulas for the intrinsic 
volumes and for the curvature measures of polyhedral cones (although the actual computation 
of the intrinsic volumes may still very well be complicated as the resulting formulas include 
volumes of spherical polytopes). 

Another class of cones, for which one has closed formulas for the intrinsic volumes, are 
smooth cones, i.e., cones C C W 1 such that the boundary M := dK of K = C PI S^ 1 is a 
smooth (i.e., C°°) hypersurface of The formulas for the intrinsic volumes involve the 

principal curvatures of M, which we shall describe next. 

In general, let M C S^ -1 be a smooth submanifold of the unit sphere. For peMwe denote 
the tangent space of M in p by T p M, and we denote its orthogonal complement in TpS^ 1 = p 1 - 
by T^-M. Let ( £ T p M be a tangent vector, and rj G T p M a normal vector. It can be shown 
that if c: M. — > M is a curve with c(0) = p and c(0) = C, and if it; : M — >■ M rf is a normal extension 
of n along c, i.e., w(t) G T^^M and w(0) = 77, then the orthogonal projection of u/(0) onto T p M 
neither depends on the choice of the curve c nor on the choice of the normal extension w of 77 
(cf. for example [42, Ch. 14] for the hypersurface case, or [14, Ch. 6] for general Riemannian 
manifolds). It therefore makes sense to define the map 

W PtV : T p M — > T p M , -U TpM (w(0)) , 

where w : R ->■ R d is a normal extension of rj along a curve c : M — > M which satisfies c(0) = p 
and c(0) = (, and IIj^m denotes the orthogonal projection onto the tangent space T p M. This 
map is called the Weingarten map. 

It can be shown that W P)J1 is a symmetric linear map (cf. [14, Ch. 6]), so that it has m := 
dimM real eigenvalues Ki(p, 77), . . . , K m (p, 77), which are called the principal curvatures of M 
at p in direction 77. The corresponding eigenvectors are called principal directions. Furthermore, 
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we denote the elementary symmetric functions in the principal curvatures by 

CT i(P> V) ■= K Ji (P> V)--- (p, v) ■ (53) 

l<jl<-..<ji<m 

When we are working with orientable hypersurfaces, i.e., with submanifolds of codimen- 
sion 1, which are endowed with a global unit normal vector field v. M — > T^M, v(p) 6 T p M, 
II^Cp)!! = 1) then we abbreviate (Ji{p) := o~i(p, f{p))- When M = dK is the boundary of a spher- 
ically convex set, and additionally a smooth hypersurface of S d ~ x , then we always consider M 
to be endowed with the unit normal field pointing inwards the set K (this implies Ki(p) > 
for all i = 1, . . . , d — 2). 

In the context of (spherically) convex sets, Weyl's classical tube formula [44] says the fol- 
lowing: Let C C R d be a closed convex cone such that the boundary M = dK of K = CD S^ 1 
is a smooth hypersurface of S^ 1 . Then, for 1 < j < d — 1, 

Vj(C) = - - • J a d ^(p) dM , (54) 

peM 

where Od-i '■= vol^-i S^ 1 = ' anc ^ denotes the volume element induced from the 

Riemannian metric on M. 

The problem is that the cones Cp >n , whose intrinsic volumes we want to compute, are neither 
polyhedral nor smooth (for n > 3). But the rank decomposition (8) yields a decomposition 
of Cp tJl into smooth pieces, which is the basic idea behind the proof of Theorem 3.5. In the re- 
mainder of this section we define the notion of a stratifiable convex set, which is a generalization 
of both polyhedral and smooth convex sets, and we state a suitable generalization of (54). 

In the following let M C S^" 1 be a smooth submanifold of the unit sphere. We may consider 
the tangent resp. normal bundle of M (cf. [40, Ch. 3]) as submanifolds of R d x R d via 

TM= [j {p} x T P M , T L M = [j {p} x T^M . 

p£M p£M 

Furthermore, if M C S" i ~ 1 , we also consider the spherical normal bundle 

T S M := \J {p} x TpM , T$ M := T^M n S^ 1 . (55) 

The tangent and the normal bundle are both so-called vector bundles, as all fibers of the 
canonical projection maps (x, v) t- > x are vector spaces. The spherical normal bundle is a sphere 
bundle, as all fibers are subspheres of the unit sphere. For the generalization of Weyl's tube 
formula we need to consider another class of fiber bundles, where each fiber is given by (the 
relative interior of) a spherically convex set. 

Let C C R d be a closed convex cone, and let He denote the canonical projection onto C. 
For p G C we define the normal cone of C in p by 

N P (C) :={veR d \U c (v + p)=p} , 

which is easily seen to be a closed convex cone with N p (C) C p 1 - . For a subset M C C, we 
define the spherical duality bundle via 

N S M := [j {p} x N^M , N^M := relint(JV p (C)) n S^ 1 . (56) 
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Note that we have not imposed any smoothness assumption yet, but ifMCCnS d 1 is smooth, 
then we have N S M C T S M. Note also that N S M in fact depends on M and C. 

Definition 4.1. Let C C R rf be a closed convex cone. We call the spherically convex set 

■ k 

K := C n stratifiable if it decomposes into a disjoint union K = [) i=0 Mi, such that: 

1. For all < i < k, Mi is a smooth connected submanifold of S"* -1 . 

2. For all < i < k the spherical duality bundle N s Mi is a smooth manifold. 

■ k 

If (1) and (2) are satisfied, then we call K = {J i=0 Mi a valid decomposition. Furthermore, we 
call a stratum Mj essential if dim N s Mi = d — 2, otherwise we call it negligible. 

The following theorem is the announced generalization of Weyl's tube formula (54) to strat- 
ified sets. A proof may be found in [6, §4.3]. Formulas similar to the one that we give in the 
following theorem may also be found in [3]. 

Theorem 4.2. Let C C such that K := C n S" i_1 is stratifiable and decomposes into the 

■ k 

valid decomposition K = \J i=0 Mi, with M$ = int(-fT) ; and Mi, . . . , M\. denoting the essential 
and Mfc + i, . . . , M~ k , k < k, denoting the negligible pieces. Then, for 1 < j < d — 1, 

k 

Vj(C) = J^CCMO, (57) 
i=i 

*i(C, Mj) = ^— -L— - • | | a^l^ip, - V ) dN*(C) dMi,fori = l,...k, 

(58) 

where di := dimMj + 2, ajp(p,—r]) denotes the Ith elementary symmetric function in the 
principal curvatures of Mi at p in direction —r\, and o~e(p, —rf) := if i < 0. 

A paraphrase of (57) is that the curvature measures vanish at the negligible pieces. Note 
also that in the sum (57) only those strata Mj contribute, for which j < di — 1. 

4.2 Orthogonal, unitary, (compact) symplectic group 

In this section we discuss the compact Lie group, which preserves the canonical scalar product 
on Fg given by (x, y) = x^y = YH=i %iVii x i V e ^ We denote this group by 

G(n) := Gf}(n) := {U G ¥f n \ Vx, y G F£ : (Ux, Uy) = (x, y)} 
= {U € F" x " | = /„} . 

For /3 = 1,2,4 the group Gp{ri) is called the orthogonal group, unitary group, and (compact) 
symplectic group, respectively, cf. for example [19, §7.2]. Note that an element U G Gp(n) may 
be identified with an orthonormal basis of by interpreting the matrix U as the n-tuple of its 
columns. We drop the index j3 to simplify the notation. 
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The Lie algebra of G(n), which we identify with the tangent space of G(n) at the identity 
matrix I n , is given by the real vector space of skew-Hermitian matrices 

Skew n := Skew^ := T In G(n) = {A G ¥f n \ = -A} . 

To specify a left-invariant Riemannian metric on G{n) it suffices to declare an M-basis of the Lie 
algebra Skew„ to be orthonormal (and then extend the metric to G(n) by pushing it forward via 
the left-multiplication). For (3 = 4 we declare the following basis of Skew n to be orthonormal: 

{lEu I 1 < i < n,i G {i,j,k}} (59) 
U {Eij -E ji \\<j<i<n}V} {i(Eij + Eji) | 1 < j < i < n,t G {i, j, k}} , 

and for j3 = 1,2 we use its intersections with M nxn and C nxn , respectively. It is readily 
checked that this yields a bi-invariant metric on G(n) (the bi-invariance in fact determines the 
Riemannian metric up to scaling). 

Before we calculate the volume of G[n) with respect to the volume determined by the above 
Riemannian metric, we recall the well-known (smooth) coarea formula. The coarea formula is 
essential for the computation in Section 4.3. 

If L : V — > W is a surjective linear operator between euclidean vector spaces V and W, then 
we define the normal determinant of L as 

ndet(L) := | det(L| ker(i) ±)| , 

where £|kcr(L) x : ker(L) -1 - — > W denotes the restriction of L to the orthogonal complement of 
the kernel of L. Obviously, if L is a bijective linear operator, then ndet(L) = | det(L)|, so the 
normal determinant provides a natural generalization of the absolute value of the determinant. 

Lemma 4.3. Let ip: M.\ — > M.2 be a smooth surjective map between Riemannian mani- 
folds Mi, M 2 . Then for any f : M.\ — > M. that is integrable w.r.t. dM.\ we have 

J fdM >= I I ^M^ d ^ { " }dM '- < 6o) 

Mi qeM2 pe<p- 1 (g) 

The inner integral in (60) over the fiber (p^ 1 (q) is well-defined for almost all q G M-2- This 
follows from Sard's lemma (cf. for example [39, Thm. 3-14]), which implies that almost all 
q G A^2 are regular values, i.e., the derivative D p ip has full rank for all p G tp~ l {q). The 
fibers (p~ 1 (q) of regular values q are smooth submanifolds of Ai\ and therefore the integral over 
tp~ l {q) is well-defined. One calls ndet(D p ip) the Normal Jacobian of at p. 

See [30, 3.8] or [17, 3.2.11] for proofs of the coarea formula with .Mi, M2 being submanifolds 
of euclidean space. See [27, Appendix] for a proof of the coarea formula in the above stated 
form. 

To get the volume of G(n) we consider its action on F§. Note that we have an inner product 
on given by (x, y) = x^y. Considering as a (/?n)-dimensional real vector space, we have 
the inner product on given by (x, j/)r := K((x,y)), where Sft: F^ — > M denotes the canonical 
projection. For (3 = 4 this amounts to the same as declaring the set {1 ■ ej | 1 G {1, i, j, k}} an 
orthonormal basis of the real (4n)-dimensional vector space H n , where G W 1 denotes the ith 
canonical basis vector. 



25 



We now get the volume of G(n) in the following way: The map ip: G(n) — > S(Fp) = 
{x G Fjg | ||x|| = 1}, U i — ^ U m 61, is a Riemannian submersion, i.e., for every £7 G G(n) the 
restriction of -Df/V to the orthogonal complement of its kernel is an isometry. In particular, the 
Normal Jacobian of 99 is everywhere equal to 1. Furthermore, each fiber ip~ 1 (x) is isometric 
to G(n — 1) as is easily checked. An application of the coarea formula thus yields vol G(n) = 
volS{¥^) ■ vol G(n - 1). As vol,S(F™) = volS^ 1 ' 1 = O fin -i, we obtain by induction 

vol <?(„) = n o«-. = n Ym = r • * ni " +mi ii, ' (61) 

i=l i=l v 2 / i=l 1 v 2 / 

In Section 4.3 we will need to consider certain subgroups of G(n). By a distribution of 
r G Z>o we understand a tuple p = (pi, . . . , p m ) G Z™ such that |p| := p± + . . . + p m = r. For 
such p with |p| < n we define the closed subgroup G(n,p) of G(n) consisting of the matrices 
having a block-diagonal form prescribed by p: 

G(n,p):={di ag (U 1 ,...,U m ,U')\U i eG(p i ),U' eG(n-r)} . (62) 

Note that G(n, p) with its induced Riemannian metric is isometric to the direct product G(pi) x 
. . . x G(p m ) x G(n—r). Furthermore, the homogeneous space G(n)/G(n, p) is a smooth manifold 
of dimension 

m m 

dim G(n)/G(n, p) = dim G(n) - dim G(n - r) - ^ dim G{pi) = dp, n - dp, n - r - ^ dp tPi 



i=i 




n — r 
2 




(63) 



The case p = \^ = (1, . . . , 1) (r-times) will be of particular importance. Note that G(l) = 
S(¥ p ) = {a G Fp I ||a|| = 1}, so that G(n, l (r) ) = Sflfy) x . . . x S'(F^) x G(n - r). We use the 
notation 

G n , r := G(n)/G(n, l< r >) . (64) 

Furthermore, we denote by G(n) — >■ G„ jr , f/ 1— >■ [{/]:= U ■ G(n, 1^) the canonical projection. 

Note that G(n) has a natural action on G n ^ r given by (Ui, [U2]) >->■ [C7iC/ 2 ] for U\,Ui G G(n). 
Moreover, as G(n) acts transitively on G„ ir , there exists up to scaling at most one Riemannian 
metric on G n<r , which is G(n)-invariant. In the following paragraphs we will give a concrete 
description of the tangent space T[ In ]G n>r , and specify on it a G(n)-invariant Riemannian metric. 

As G(n,lW) = I (0 u') A = dia g( A i'---' A r)> A i e ^(F^), t/ r G G(n-r)|, the tangent 
space of G n!r at [In] is given by 

T In G(n, lW) = I ^ 5) = diag(ai, . . . , a r ) , »(oi) = , fit = _ 5 l . 

The orthogonal complement of T/ n G(n, l( r )) in Tj n G(n) = Skew n , the space of skew-Hermitian 
matrices, is given by 

m^:=(T In G(n,l {r) )) ± = l [ (y | X G *t = -X, Y G F^*'} . (65) 
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It can be shown (cf. [25, Lemma II. 4.1]) that there exists an open ball B around the origin 
in Ti n G(n) = Skew n such that the intersection I?nSkew n is diffeomorphic to an open neighbor- 
hood of [I n ] in Gn )T . Moreover, the tangent space of G n)T in [I n ] may be identified with Skew n , 
and the restriction of the inner product on Skew n to Skew n yields a well-defined Riemannian 
metric on G n ,r, which is G(n)-invariant. (See [7, §5.2] for a more detailed description of the 
induced Riemannian metric on a homogeneous space in a similar situation.) 

For j3 = 4 we have the following orthonormal basis of Skew n , cf. (59), 

{Eij - Eji | G X} U {i(Eij + Eji) | i G {i,j,k}, G X} , 
where X := X\ U X2 with 

Xi :={(i,j) \ l<j<i<r}, X 2 := {(», j) \ r + 1 < i < n , 1 < j < r} . (66) 



For further use in Section 4.3, we denote this orthonormal basis of Skew n = T[ In ]G n , r (for 
P = 4) by 

Vij ■= E ij ~ E H , Vij ■= i>(Eij + Eji) , G X, t G {i, j, k} . (67) 

For (3 = 2 we have the orthonormal basis {77^ | 1 G {1, i}, G X}, and for f3 = 1 we have the 
orthonormal basis {77^ | (i,j) G X}. 

The canonical projection G(n) — > G n)T turns out to be a Riemannian submersion, and an 
application of the coarea formula yields 

volG(n) volG(n) = volG(n) 

n ' r volG(n,lM) {vol S(¥ P )Y- vol G(n-r) Oy x ■ volG(n - r) ' 1 J 

4.3 Deducing the formulas for &j(/3,n,r) 

In this section we derive the claimed formulas for the curvature measures $j(/3,n, r) stated 
in Theorem 3.5. The organization is straightforward and is as follows: First, we recall the 
well-known face structure of Cp in , which is described for example in [8, §11.12]. Then we show 
that Cp jU is a stratified cone, and we determine the essential and the negligible pieces. Finally, 
we compute the principal curvatures of the strata, so that we can give the proof of Theorem 3.5 
by an application of Theorem 4.2 and a small computation. 

In this section we change to the spherical viewpoint and denote the intersection of Cp >n with 
the unit sphere by 

K n := C^ n n S(Rer^ n ) = {A G Her^ \AhO, \\M = 1} . 

First, we recall the well-known face structure of Cg >n , cf. Section 2.1. 

Proposition 4.4. The faces ofCp^ n are parametrized by the subspaces o/F^. More precisely, 
for LCI^ an F 'p -linear subspace of ¥p -dimension r, where < r < n — 1 the set {A G Cp jTl \ 
im(A) C L} is a face of 'Cp <n of dimension dp >r . Conversely, every face ofCp^ n is of this form 
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and hence has dimension dp )T for some < r < n — 1. Moreover, every face of Cp jn is of the 
form 

U { A o) Uj A ' eC ^} , (69) 
where U 6 Gp(n) and < r < n — 1. The normal cone at the face defined in (69) is given by 

(70) 



Mo -A>) U '\ A " eC ^- r } ■ 



Proof. See for example [8, §11.12]. The proof for the real case given there extends to the complex 
and the quaternion case in a straightforward way. □ 

Note that (69) and (70) show that when analyzing a face oiCp tTl , by choosing an appropriate 
basis of F§, we may assume without loss of generality that this face is of the form Cp >r x {0} 
with corresponding normal cone {0} x (—Cp, n - r ). 

Next we will show that K n is stratifiable and exhibit a valid decomposition of K n , cf. Def- 
inition 4.1. For this, we use a finer distinction than in (8): we classify the matrices A e K n 
not only according to their ranks, i.e., according to the multiplicity of the eigenvalue 0, but 
also according to the multiplicities of the nonzero eigenvalues. To achieve this we define the 
eigenvalue pattern of an element A G K n via 

patt(A) := (pi, . . . , p m ) , iff Ai = . . . = A Pl > A pi+ i = . . . = A Pl+P2 > . . . , 

where Ai > . . . > A r > are the positive eigenvalues of A. Note that patt(A) is a distribution 
of r = ik(A). The spherical cap K n thus decomposes into 

■ n 

K n = LLiUl \= r Mn ^ ' Mn >? := {A G Kn 1 P att (^) = P} ■ ( 71 ) 



n-1. 



Note that mt(K n ) = \J\ P \= n M n, P and dK n = (J r=1 U\p\= r M n, P - 

Proposition 4.5. T/ie set M HjP , \p\ < n, defined in (71) is a smooth submanifold of the unit 
sphere S^Her^). Furthermore, the duality bundle N s M n ^ p defined in (56) is a smooth manifold 
for all \p\ < n. Hence (71) is a valid decomposition. 

The strata {M lW | 1 < r < n) , where l( r ) := (1, 1, ... , 1), are essential and all the other 
strata M„ iP are negligible. 



Proof. We fix a pattern p = (pi, ... , p m ) with |p| = r < n, and define the set P C S n via 



A e S 71 ' 1 



Ai — . . . — A Pl > A Pl+ i — ... — A Pl+P2 > . . . 

. . . A^. ^> \r-\-l — ••• — A77, = 



(72) 



It is easily seen that the set P is the intersection of the interior of an m-dimensional polyhedral 
cone with the unit sphere, so it is an (m — l)-dimensional submanifold of S n ~ l . We now consider 
the map 

Y> : P x G{n) -> Her^ in , (A, U) ^ U • diag(A) • f7 t . (73) 

This map is smooth, and its image is M n ^ p by the principle axis theorem. Concerning the fiber 
of A G M n ^ p , we may assume w.l.o.g. that A = diag(p). Note that U ■ diag(A) • = diag(p) iff 
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X = p and U • diag(A) = diag(A) • U. Furthermore, it is easily checked that for A G P, we have 
U ■ diag(A) = diag(A) ■ U iff U G G(n,p). We may thus conclude that, for any A G M raiP , 

V>(A, 17) = A ^ A = // and U G G(n, p) . (74) 

This implies that the map tp factors over the product P x (G(n)/G(n, p)), i.e., we have a 
commutative diagram 

P x G(n) > Her^ 



ii 

Px(G(n)/G(n,p)) 



1 • < 75 > 



where II: PxG(n) — > Px (G{ri)/G{n, p)) denotes the canonical projection map. Moreover, (74) 
implies that the map V is injective, and thus a bijection on the image of tp, which is M njP . It 
is straightforward to show that the derivative of ip has full rank. Since P x (G{n)/G(n, p)) 
is compact, it follows that tp is a homeomorphism onto its image M„ ]P . Hence M„ iP is a 
submanifold of Herg„ and ip induces a diffeomorphism of P x (G(ri)/G{n,p)) to M njP , cf. for 
example [10]. 

As for the claim about the duality bundle, note that for A = U ■ diag(A) ■ U\ A G P, the 
normal cone NA(K n ) is given by (cf. Proposition 4.4) 



^4 G Cp >n — r 



Denoting by int(TT n _ r ) the interior of ifn-r with respect to the topology on S^Her^-j.), we 
can define the map 

*: P x G(n) x int(7C n _ r ) ->• Her^ x Her^ , 

(A,C/,A")^ (^(A,C/), u(l _V)^ f ) . 

This map is smooth, and its image is N s M n)P , cf. (56) and (70). If *(A, U,A") = *(p, V,P") 
then ip(X, U) = ip(p, V), i.e., A = p and V = U • D for some D G G(n, p). Furthermore, for the 

second component we get U ^ = UD b'^J D^U^, i.e., 

A" = D' • B" ■ (D')t , (76) 

where D = diag(Z?i, . . . ,D m ,D'), cf. (62). Using the conjugation (76) to define an action of 
G(n, p) on mt(K n - r ), we may form the factor space (G{n) x int(7C n _ r )) /G(n, p). As the above 
arguments about the injectivity of ^ are reversible, it follows that * factors over the product 
P x (G(n) x mt(K n - r ))/G(n, p). We thus obtain the smooth injective map 

P x (G(n) x mt(K n - r )) /G(n, p) -»• Her^ x Her^ , 
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whose image is given by N s M n>p . Arguing as above shows that N s M n>p is indeed a smooth 
manifold. Its dimension is given by 

dim N s M n)P = dim P + dim G(n) + dim K n ^ r - dim G(n, p) 

m 

= (m - 1) + dp >n + {dfa n - r - 1) - (dp, n -r + ^ d P'Pi) 



i=l 



Note that since m < r, we have 



m / \ 



= if m = r, i.e., p = l( r ) 
< otherwise . 



Therefore, as dimHerg n = dp >n = n + /3 • (™), we have 



dim N b M n , p 



= dimHer^-2 ifp = l( r ) 
< dim Her ^ — 2 otherwise. 



Therefore, the strata M n 1 ( T .), 1 < r < n, are essential, and all other strata are negligible. □ 

From now on we may restrict the computations to the essential strata M n ^ r ). We use the 
notation 

P r := {A e S' 1 ' 1 | Ai > A 2 > . . . > A r > = A r+i = . . . = A n } . (77) 
From the proof of Proposition 4.5 we obtain the following corollary. 

Corollary 4.6. Let 1 < r < n, and let <p r be defined via 

<p r :P r x G n , r -> M n>lW , (A, [U]) ^ U ■ diag(A) • C/t . (78) 
Then ip r is well-defined and a diffeomorphism. □ 

In the following proposition we compute the Normal Jacobian of (p r , i.e., the absolute value 
of its Jacobi determinant. This result is well-known, but we include the proof for the sake of 
completeness. 

Proposition 4.7. The Normal Jacobian of <p r at (A, [U]) is given by 

ndet(D (A>|p]) = r<*~-W* . A(A) /3 • f[ X^ r) , 
where A(A) = Y\i<i<j< r i^i ~ -\?) denotes the Vandermonde determinant. 
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Proof. We restrict to the quaternion case f3 = 4; the cases (3 G {1,2} are similar. By orthogonal 
invariance, we may assume w.l.o.g. that U = I n . For ( G T\P r we have 

I>(A 1 [/„])yr(C,0)=diag(C). (79) 

Recall that we identify T^ n ]G n>r with Skew n , which has the orthonormal basis (rj^) defined 
in (67). Let U-j : R — > G(n), with t G {l,i,j,k} and (i,j) G X (cf. (66)), be curves such that 
the induced curves [U-j] : M — >■ G njr define the directions j^-, cf. (67). Then we may compute 
the derivative of (p r in the second component for i = 1 by 

£W„])¥>r(0,»&) = |(aj(t) • diag(A) • C/J(t)t)(0) 
= ^ • diag(A) - diag(A) • 77^. 



(80) 



(Aj - Ai) • (^ + Eji) if 1 < j < i < r 



[Xj ■ (Eij + E^) \ir + l<i<n, 1 < j <r , 

and similarly for 1 G {i, j, k} 

£>(A,[/„])¥'r(0,^) = fM 3 {t) ■ diag(A) • C^-(t)t)(0) 
= Vij ■ diag(A) - diag(A) • r)\j 

= U\j - Xi) ■ L(Eij - E^) if 1 < i < i < r 
" |Aj ■ t^jj - Eji) ifr + l<i<n,l<j<r. 

Note that the D^ x ^ In ^ip r (0,r]^) are pairwise orthogonal elements of Her^. Moreover, the 
directions VijiVljirfijiVij have length 1, while the directions E^ + Eji,i(Eij - Eji),j(Eij - 
Eji),h(Eij — Eji) have lengths y/2. Taking this into account, we get 

ndet(D (A)[u])¥ v) = II (V2-(A,-A,)) 4 - ]J (y/2 ■ A,) 4 



l<j'<j<r 7'+l<i<n 

i<j'<r 



= 2 r ( r - 1 )- ]f (A, - A,) 4 • 2 2r ^ ■ H Xf n ' r) . □ 

l<j<i<r 

It remains to compute the principal curvatures of the essential strata M ^ before we can 
use Theorem 4.2 to compute the intrinsic volumes of K n . We will do this in the following 
proposition. Note that the dimension of M is given by 

dimM n lW = dimP r + dimG n , r { = ] r - 1 + (3 • ((") - ("7)) 

= Pr(n-r) + r-l + P($ ■ 

Proposition 4.8. Let A = [/-diag(A) • W G M n 1 ( r ) with X G P r . Furthermore, let A" G Cp, n - r , 

so that B := U ■ ^ • C/^ G Na{K u ) is a vector in the normal cone of K n at A. If 

fii > . . . > fi n - r > denote the eigenvalues of A" , then the principal curvatures of M n 1 ( r ) at A 
in direction —B are given by 

— — — , — - — , • • • , — - — (each of these values p -times J , 
Ai Ai A2 A2 A r A r 

0, . . . , (overall (r — 1 + /3 (Q) -times) . 
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Proof. By orthogonal invariance we may assume w.l.o.g. that U = I n , so that A = diag(A) and 
A" = diag(/x). From Proposition 4.7 we get that the tangent space of M n 1 ( r ) at A is given by 
(omitting the argument (A, [I n ])) 



T A M nlW = D<p r (T x P r x Skew n ) . 

It is easily seen that all the vectors in Dip r (T\P r x {0}) are principal directions of M n 1 ( r ) at A 
with principal curvature 0, thus giving r — 1 of the claimed r — 1 + zero curvatures. 

Concerning the second component, we again only consider the quaternion case /3 = 4, 
the other cases being similar. As in the proof of Proposition 4.7 let U-j : R — > G(n), with 
i G {1, i, j, k} and G X (cf. (66)), be curves such that the induced curves [U^] : R — > G n>r 
define the directions 77^, cf. (67). We denote the images of Dip r by 

Qj -.= D<p r (o, Vij) e r A M n>1(r) , «, g {1, i, j, k} . (82) 

Explicit formulas for the vectors Qj are given in (80). 

We define normal extensions of —B = diag(0, jj) along the curves ip r (\, [U-j(t)]) via 

v\j(t) := E^-(t) • diag(0,^) • U^tf , t G {l,i,j,k} . 
Differentiating these normal extensions t = yields for i = 1, using rj^- = E^ — Eji, 

|4(0) = {Eij - Eji) • diag(0, n) - diag(0, M ) ■ (E tj - Eji) 
Jo ifl<j<i<r 
I -fH-r • (-Eij + Eji) if r + 1 < i < n, I < j < r , 

and for 1 G {i, j, k}, using r?^- = t(£' ij - + E^), 

|^(0) = i{Eij + E^-j) • diag(0, /i) - diag(0, M ) ■ + 

Jo ifl<j<i<r 

I -/Xi_ r • i(Eij - Eji) if r + 1 < i < n, 1 < j < r . 

Comparing this with the values of (ij,Qj,(jj, Ctj gi ven i n (80), implies for 1 G {1, i, j, k} 

-„*.((,) = J°'^ ifl<i<^<r 

dt 13 y ' \- >J tf L -Qj if r + 1 < i < n, 1 < j < r . 

We conclude that the directions Qj , Qj , ■ , C^. are principal directions with curvature and , 
respectively. □ 

Before we give the proof of Theorem 3.5, we state a small lemma, that will come in handy 
for an integral conversion later on. 

Lemma 4.9. Let f : W l \ {0} — > R be a homogeneous function of degree k, i.e., f(x) = \\x\\ k ■ 
/(||x|| _1 ■ x). Then for a Borel set U C S^ 1 and U = {s ■ p | s > , p G U} 

f If \M 2 

J f(p) d p = 2 ^ ^ • J e_ 2 -mdx. 

pel/ { 2 ' x& (j 
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Proof. The normal Jacobian of the projection II: W l \ {0} — > S n 1 , II(x) = ||x|| 1 • x, is given 
by ndet(D x II) = ||x||-( n_1 ). From the coarea formula (Lemma 4.3) we thus get 



e a 



f(x)dx = J J s n ~ l -e-" 2 / 2 ■ f(s-p)dsdp= J f(p)dp- J s n ~ 1+k -e- s2 / 2 ds. 



I , , , 

J J JO J JO 

xeu v& v& 

Substituting t := s 2 /2, and using the well-known formula J °° t z ~ l ■ e~* dt = Y[z), we obtain 

poo , /*oc 

/ S n - 1+k .e- s2 ' 2 ds = 2^- 1 . t n ¥- 1 -e- t dt = 2 n i k - 1 .T(^). □ 
Jo 

Proof of Theorem 3.5. Recall from (9) that $>j(/3,n,r) = ®j(Cp jn ,W/3 jnjr ) and in the stratifi- 
cation (71) of K n = Cp :n n S^Her^) only the strata M n 1 ( r ) are essential, cf. Proposition 4.5. 
Denoting <5(n, r) := dimM n 1 ( r ) + 2 = j3r(n — r) + dp tT + 1, we thus obtain from Theorem 4.2 

= p. .i / / ^i,r)-i-l(^- S ) ^ dM n,l« > 

AeM n,iW Se7V| 

where the superscript in cr^ shall indicate the dependence on M n 

Recall that in Proposition 4.8 we have computed the principal curvatures of M n 1 ( r ). Using 
the notation (x y _1 ) x/3 = (x . . . , x (/3-times) and x y _1 defined in (38), we 
obtain 

*j(P,n,r) = - —± — j J a (5(n , r) _ j _ 1 ((A" 1 M ) x/3 ) dtfj dM njlW , (83) 



where A and // denote the (positive) eigenvalues of ^4 and — B, respectively. Using the relation 
• • • ' ^) = ( Xl ' -- x n)~ N - °N-k(xi, ■ ■ -,xn) and observing 5(n,r) - j - 1 = /3r(n - r) + 
rf/3,r — ji we can rewrite the integrand via 



lli=i Aj 

In Proposition 4.7 we have shown that the Normal Jacobian of at (fx, [U]) is given by 
2«(n-i)/3/4 . A(^i) /3 . The normal cone of at ^4 G M n l(r) is isometric to Cp, n - r , cf. Propo- 
sition 4.4. Further, M n _ r y( n -r) equals K n _ r = Cp >n - r D S^Her^^) up to strata of lower 
dimension. Taking into account the Normal Jacobian of <£> n - r determined in Proposition 4.7 
(note that n needs to be replaced by n — r), we can transform the inner integral of (83) via the 
coarea formula to obtain 

j f(\^)dNi= j /(A,m)-2(-^-'- 1 ^/ 4 -A( / x)^( M ,[C/ 2 ]), 

BeJVf P„- r xG„- 



-r.n — r 



where we have abbreviated /(A, /x) for the integrand (84). 
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Similarly, we may transform the outer integral of (83) by applying the coarea formula to 
the map <p r . As a result we obtain 

j j f(\,ridN%dM nil{r) = j j f{\,n) 

AeM n,iM B^N% (A,[l/i])eP r xG„, r (jl, [U 2 ] ) G P n -r x G n _ r ,„_ r 

• y^-r-W* • A(Xf • fl Af ( "- r) • 2^- r ^ n - r - 1 ^ 4 • A^f d(fi, [U 2 ]) d(X, [Ui]) . 

1=1 

(85) 

Note that we have 

(68) volG(n) vol G(n — r) vol G(n) 



vol G n ,r ■ vol G n - 



r,n—r 



Oy x • vol G(n-r) 0™Z[ 0^ ' 

Replacing f(X,fi) again by (84), the integral (85) simplifies to 

ft 1 p p 1 

-* r r n — r 

Fn(/3 / 2 ) ' ' • y y A ^ • A ^ • fpj-bA**) dxd ^ . ( g6 ) 

-*r r n—r 

where in (*) we have used (37) and the small computation 

VQlG(n) (61) 2" • W"+1W4 ■ n? =1 ^ = n(n _ i)/?/4 » (4J) ff n(n-lW4 , „, 

^-i (2-^/2/r(f))« " fir(f) F n (y3/2) 

The integrand in (86) is bihomogeneous in A and /j,. Its degree in A equals fify +j — dp jr = j — r, 
and its degree in \i is given by Pi^ 1 ^) + /3r(n — r) — j + dp^ r = — j + r. Using Lemma 4.9 
twice, we get 

, oc . (2 7 r)™(«- 1 W 4 • n! 2 2 

(86) - - 



K(/?/2) 2J/2.r(|) 2 (p(l)- j+ n)/2 . r ^ 



+n-i > 



/ y e-^^^ • A(A) /3 • A( M )" • U d - dptr (X;^) dXdfi 



Pr Pji — r 



F n ((3/2) ■ (2tt)»/2 

P P 



/ / e - M ^-A(A)^A(^./,,_^ r (A; / ,)cZA^ 



r J n — r 



The positive orthant decomposes into r! isometric copies of P r , such that their interiors are 
disjoint. More precisely, the copies of P r are parametrized by the permutations of {1, ... ,r}, 
which indicate the order of the components of a vector in K_. The same applies to M. r t~ r 



34 



and P n -r- As the Vandemonde determinant is antisymmetric, and //3,fc(A; //) is symmetric both 
in A and in fi, we finally obtain 

*Mn,r) = E ■ J J e- M ^-A(A)^A(^. //3j _^(A; M )dA^ 



Pr Pfi — r 

-(:) >, W! |'.(2,)^ / e -*-|A(A)P-|A(,)r/ W -^(A;,) rf , 

i/:=(Ajt)eK!£ 

(40) AA J/3{n,r,j - dp, r ) n 
W ' " F n ((3/2) " • 

A Support measures 

The kinematic formula, as we state it in (22) and which we need in the proof of Theorem 1.3, 
follows, as well as the other formulas in Theorem 2.3, from a general kinematic formula for 
support measures that Glasauer has proved in [20], cf. also [21, 28]. As this general formula is not 
easy to trace in the literature, we will give in this section a short account of support measures, 
and we will derive Theorem 2.3 from Glasauer's kinematic formula for support measures. 



A.l Definition and properties 

A major drawback of the curvature measures is that while they localize the notion of intrinsic 
volumes, they do not possess a duality structure. In fact, one may interpret the curvature 
measures as primal localizations of the intrinsic volumes, while the support measures may be 
thought of as primal-dual localizations of the intrinsic volumes. They generalize the curvature 
measures, and they do admit a duality structure. To obtain this duality we need to replace the 
Borel algebra 3S{£), defined in (1) by another cr-algebra, which we define next. 

Definition A.l. Let £, £' be finite-dimensional euclidean spaces, and let 33(£ x £') denote the 
Borel cr-algebra on £ x £'. Then we call 

£') := {M G 38{£ x £') | VA, fi > : (A, y)M = M} , (87) 

where (A, fi)A4 := {(Xx,fiv) \ (x,v) G Ai,x G £,v G £'}, the biconic (Borel) a-algebra 
on (£,£'). If £' = £ then we call £$(£,£) the biconic (Borel) a-algebra on £. 

It is straightforward to show that the biconic cr-algebra 3${£,£') satisfies the axioms of 
a cr-algebra. In the following we will only consider the case £' = £. In this case we have 
the natural involution .* : £ x £ — > £ x £, (x,v)* := (v,x). This map induces an involution 
M i— > M* :={(v,x) | (x, v ) G Ai} on &(£, £). The group of isometries of £, i.e., the orthogonal 
group 0(dim£) having fixed an orthonormal basis in £, acts on &(£, £) via QM := {(Qx, Qv) \ 
(x,v) G M}. 

The support measures of C are measures on the biconic cr-algebra on £ 

e (C, .),..., e d (C,.): &(£,£) . 

In the following we will only give the definition of the support measures for polyhedral cones 
in analogy to (14) and (15). By continuity one obtains from this also a definition for support 
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measures for general closed convex cones. More precisely, it can be shown, cf. [20], that 0., 
can be extended to the set of all closed convex cones in such a way that if the sequence 
Ci, C 2 , . . . converges in the Hausdorff metric to C, then lim^oo @j(Ck, M) = @j(C, M.) for all 
M G 3S(£, £). We skip a definition via volumes of tubes similar to Definition 3.2 as we will not 
need such a characterization. 

Let C C £ be a polyhedral cone in (i-dimensional euclidean space £. It is well-known, 
cf. for example [46, §2] or [22, §14], that the set of faces of a polyhedral cone forms a lattice. 
Furthermore, the face lattice of the dual C is the opposite lattice of the primal cone C, cf. [46, 
§2.3]. More precisely, if F C C is a face of C and if F := remit (F) denotes the relative interior 
of F, then we denote the corresponding face of the dual cone C and its relative interior by 

F° := span(F) ± nC, F° : = relint(F <> ) . 

Let J-j denote the set of the relative interiors of the j-dimensional faces of C. If lie an d 14^ 
denote the canonical projections on C and on C, respectively, then the support measures of C 
evaluated in M G &(£,£) are given by 

e,(c,.M)=V Prob [(n c (x),n a (x)) e (Fx^)n^l , j = o,i,...,d. (88) 

For j > d we define &j(C, M) := 0. 

Analogous to Proposition 2.1 we formulate in the following proposition some well-known 
facts about the support measures, which may be verified using the characterization of @j in 



Proposition A. 2. 1. The support measure &j(C, .) is concentrated on C x C, that is, 
Qj(C,M) = &j(C,M n (C x C)). Furthermore, we have Qj{C,M x C) = <f>j{C,M) 
forMe &(£). 

2. The support measures are invariant under orthogonal transformations, i.e., for Q G 0(d) 
we have @j(QC, QM) = Qj(C,M). 

3. For the support measures of the dual cone we have Qj(C,M) = @d-j{C,M*). □ 
A. 2 Lattice structures 

In this section we further examine the structure of the set of closed convex cones and of the 
biconic a- algebra. This will provide a formal framework for the important duality structure, 
which is a specific property of spherical convex geometry. 

Let the set of closed convex cones in euclidean space £ be denoted by 

:= {C C £ | C closed convex cone} , 

and let the operations A, V be defined on ^(S) via 

-.C :=C , C\ A C 2 := Ci D C 2 , C x V C 2 := C x + C 2 = cone(Ci U C 2 ) , 

where cone M := {Ax + fiy \ A, /U > 0, x, y £ M} for M C 5. It is easily checked that the struc- 
ture ), A, V, £ , {0}, -i) satisfies the axioms of an orthocomplemented lattice, cf. [9], where 
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the elements £ and {0} are the neutral elements w.r.t. the operations A and V, respectively. 
In particular, we have the de Morgan's laws 

dual(Ci n C 2 ) = C\ + C 2 , dual(Ci + C 2 ) = C l nC 2 . (89) 

The projection on linear subspaces of £ fits in this framework in the following way. Let 
C G ¥?(£ ) and let W C £ be a subspace of codimension m. Then we have 

dual(C n W) = -.(-iC A W) = C V -w = C + W L n w (C) x M m . (90) 

So projection on linear subspaces is dual to intersection in a natural way. 

The biconic er-algebra &(£,£') has far less structure than c €{£). Nevertheless, we may 
define operations A and V on {%{£,£') in the following way. For M G &(£,£'), x G £ , v G £ ', 
we introduce the fibers 

M x :=ILi\x)nM , M v :=U 2 l {v)r\M , 

where IT : £ x £' — > £ and II2 : £ x £' — > £' denote the canonical projections onto the first and 
the second component, respectively. We define the operations A and V on £$(£,£') via 

M A A := (J ({x} x (M x +M x )) , M V A := [j ((Af +M V ) x {w}) , 
xe£ veS' 

for Af, A G £'). Note that for products M = Mi x M 2 , A = Ni x N 2 G =#(£, 5') we have 

(Mi x M 2 ) A (JVi x 7V 2 ) = (Mi n N ± ) x (M 2 + JV 2 ) , (91) 
(Mi x M 2 ) V (JVi x N 2 ) = (Mi + Ni) x (M 2 n N 2 ) . (92) 

The following identities are readily checked for M,Af,0 G &(£,£'): 

MaM = M AM , M A(A AO) = (M AAA) AO , 

Al V AT = AT V M , AlV(Av0) = (AlvA)V0. 

Furthermore, the elements £ x {0} and {0} x £ ' are neutral w.r.t. the operations A and V, i.e., 

M A (£ x {0}) = M , M V ({0} x £') = M . 

In the case £ = £' , to which we will restrict ourselves from now on, we also have an involution, 
which we shall denote by given by 

->M := AT = {(«, x) I (x, u) G Al} . 

With this involution we also have identities resembling de Morgan's laws 

(M A A)* = Al* V A* , (Al V A)* = AT A A* . (93) 

Unfortunately, £) is no£ a lattice with respect to the operations A and V, the reason being 

, s in general , s 

M AM = (J ({x} x (M X + M X )J / |J ({x} x A^J = M , 



37 



as the fiber M x need not be convex. 

For a cone C G ^(£) and a biconic set M G ^(£, 5), we abbreviate -i(C, M) := (->C, = 
(C, M*). In particular, we may rewrite the property of @j stated in Proposition A. 2(3) in the 
form 

e j (C,M) = @ d - j (^C,M)) . (94) 

Furthermore, we abbriviate (C,M) A (D,jV) := (C A D, M A A/") and (C,A4) V (D,N) := 
(CvD,MvM). 



A. 3 Glasauer's kinematic formula for support measures 

We let the orthogonal group 0(d) operate on pairs (M,J\f) of biconic sets M, N G £§(£,£) via 
Q(M,Af) := (QM,QM) for Q G O(d). The following theorem is Satz 6.1.1/6.1.2 in [20] and 
Theorem 9/10 in [21], cf. also [28, §2.4]. 

Theorem A. 3 (Glasauer). Let £ be a d- dimensional euclidean space, let C,D G ^(£), and 
let M,M G £%(£,£) be such that M C C x C and N <^ D x D. Then for uniformly random 
Q G 0(d) and 1 < j < d - 1, 



E 



d-l 

@ j {(C,M)AQ(D,Af))]= ®k(C,M) ■ @ d+j - k (D,J\f) , 

fc=j+i 

Qj((C,M) V Q(AAA))J = ^2Q k (C,M) ■ @j- k (D,J\f) . 

k=i 



(95) 



(96) 



Remark A. 4. The two formulas (95) and (96) are in fact equivalent, as one follows from the 
other via duality. For example, assuming (95), we obtain 



E 



Q j ((C,M)VQ(D,M)) 



(94) 



E 



Q d - j (^((C,M)VQ(D,Af))) 



(89)^(93) 



(94) 



Q 
d-i 



e d - j (^(C,M)AQHD,Jtf))) 



(95) 



d-l 



^ e d - k (c,M)-e k+j - d (D,X) 
k=d -i+i 



e k ^(C,M))-@2 d - j -kh(D,M)) 

k=d-j+l 

^i-^^e^c^-e^AAT). 



We finish this section by deriving Theorem 2.3 from Theorem A. 3. 

Proof of Theorem 2.3. Let C C M. d be a closed convex cone, and let W C M d be a uniformly 
random subspace of codimension m G {1, . . . , d — 1}. Then for M G ^(M rf ) with M C C and 
1 < j < d - m we have 

*j(C n vf, m n vf) Prop -> 2(1) @j (c n w, (m n vf) x duai(C n w)) 

{ = } e^c n w, (M n w) x (c + vf x )) = ) @j(c n w, (m x c) a x w^)) 
= e j ((c, M x (7) a (w, w x vf x )) . 
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Writing W = QWq for some fixed linear subspace Wo C R d of codimension m and with Q G 0(d) 
uniformly at random, we thus obtain 



E[$j(Cn W,M)] =E[6 i ((C,M x C)AQ(W ,W x wj-))] 

d-i 

(95) 



W V e k (c, Mxd)- e d+j - k (w , w x w^) = <^- +m (c, m) , 

, ^ .. ' v ' 

k=J+l 

which is (20). This implies (21) via 



k + 1 =*fe(C,M) , , , Prop. 2.1(2) p 

y ' =V d+] - k (R d - m ) P = 's k . 



d—m 



E[v (cnw)] rop =- lj E i-^^(cnw) = 1 — E e [Vj (c n w)] 

3=1 J J"=l 



d—m 



d—m 

( = 1 " E W<?) Pr ° P = ^ ^o(CT) + Vi(C7) + . . . + U m (C) . 

J'=l 

As for the projection formula, we have 

^■(n w (C),n w (M)) ^-J- 1 ^ $ J+m ( C + w ± ,M + w ± ) 

Prop.^.2(l) 0j . +m(c , + ^ (M + ^ x dual(c , + w ±^ 

{ = ] e J+m (c + w ± ,(M + w ± ) x (cnw)) 

{ = } Q j+m (C + W ± ,{M x C) V (W x , W x x W)) = @ j+ m{(C, MxC)V (w x , w x x w)) . 
It thus follows that 

E[^(U W (C),U W (M))] =E[@ j+m ((C,Mx C) VQ(W ± ,W ± x W ))] 

j+m-l 

( = } V e fc (c, m x c) • e^-fcCWo 1 , <- x w ) = ^(c, m) , 
ifci v v ' s " ' 

=**(C,Af) =y j+m _ fc (]R-)=5 fcJ 

which is (22). Analogously to the above computation we obtain (23). □ 
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